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Abstract 

It has been generally assumed, since the work of von Karman and Howarth in 1938, that 
free decay of fully-developed turbulence is self-similar. We present here a simple phenomeno- 
logical model of the decay of 3D incompressible turbulence, which predicts breakdown of 
self-similarity for low-wavenumber spectral exponents n in the range ric < n < 4, where ric is 
some threshold wavenumber. Calculations with the eddy-damped quasi-normal Markovian 
approximation give the value as ric ~ 3.45. The energy spectrum for this range of exponents 
develops two length-scales, separating three distinct wavenumber ranges. 
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The decay of homogeneous, isotropic turbulence is a classical problem. For a good review, 
see Lesieur^ (chapter VII, section 10). Since the early work of von Karman and Howarth^, it has 
often been assumed that the decay is self-similar for length-scales outside the dissipation-range. 
Precisely, this assumption means that the time-dependent energy spectrum may be written as 

E{k,t) =v'^{t)£{t)F{ke{t)), (1) 

where v{t) is the rms velocity fluctuation, £{t) is the integral length-scale, and F{k) is a di- 
mensionless scaling function (see section 10.2 of Lesieur^). However, recent studies of two 
exactly soluble models — the Burgers equation^ and the Kraichnan white-noise passive scalar 
equation^ — have shown that such self-similarity does not always hold. In particular, Gurbatov 
et al.^ observed that decaying Burgers turbulence develops two distinct length-scales when the 
low wave number spectral exponent n lies in the range 1 < n < 2. The energy spectrum can 
then no longer be divided into just a low-wavenumber range k£{t) ^ 1 with E{k,t) ~ Ak'^ 
and an inertial range ki{t) 3> 1 with E{k,t) k~^. Instead a new spectral range develops 
intermediate to these two with E{k, t) ~ C{t)P, C{t) oc (i — to)^^^ ln~^/^(t - to). Since it was 
the first author of reference 3 who observed this state of affairs in Burgers decay^, we call it 
the "Gurbatov phenomenon". 

The explanation of this new range lies in the phenomenon of a k^ backtransfer for Burgers 
dynamics, the analogue of the A;^ backtransfer discovered by Proudman and Reid^ in 1954 
for the 3D Navier-Stokes equations. According to traditional beliefs, the backtransfer term in 
Burgers should be overwhelmed at low-wavenumbers by the original Ak^ spectrum, which, for 
1 < n < 2, is asymptotically much the larger. This statement, however, ignores the fact that 
the coefHcient C{t) of the backtransfer term is growing in time, while the coefficient A of the 
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lowest wavenumber spectrum ~ fc" is independent of time. This leads in the Burgers decay 
to a new length-scale i^{t) ^ with the k'^ backtransfer spectrum dominating throughout 
the intermediate range M.^{t) ^ 1 and kl{t) <^ 1. An analogous "Gurbatov phenomenon" was 
found to occur in the Kraichnan passive scalar model^. It is our purpose to present a similar 
theory of the breakdown of self-similarity for the 3D incompressible Navier-Stokes equations. 

Let us consider first the case in which the energy spectrum is dominated by the A;^ back- 
transfer term at the lowest wavenumbers (as it will be if the initial spectrum has n > 4). In 
this case, one may suppose that 

{C{t)k^ ki{t) < 1 

(2) 
e2/3(t)A;-5/3 ki{t) > 1. 

The continuity of the spectrum at the juncture k£{t) = 1 imposes a relation 

e{t) oc C^/^{t)r'''/\t). (3) 
Two other relations follow by standard Kolmogorov dimensional analysis: 

and 

Needless to say, all three of these relations hold in the standard analytical closures such as 
the eddy-damped quasi- normal Markovian approximation (EDQNM). Because we have three 
relations and three unknowns {C,i,e), we can find a solution. In fact, using (^) to eliminate 
e(t), we obtain from the other two equations 

^(t) ^ cyHt)ry\t) (4) 



and 




(5) 



We can then derive from (^) and (^) that 



dC C 




so that, for some arbitrary power p, 



Cit)cx£P{t). 



(6) 



Unfortunately, it does not seem possible to calculate the precise value of p without making 
explicit use of the Navier-Stokes dynamics, either exactly or within a closure approximation. 
We shall quote below the value which follows from EDQNM results. We can only say a priori 
that we expect p > 0, since C{t) should grow with time because of the backtransfer. In what 
follows we shall not specify p, but work with an arbitrary value. 
Eliminating C{t) by substituting (|6|) in (^) gives 



-(i)o<f-^(*) 



and hence 



lit) (X{t- to) 7-P 



(7) 



and 



2p_ 



C{t) (x{t- to)^-p. 



Notice C{t) oc £'P{t) substituted into gives also 



e{t) oc £^~{t). 
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The law of the energy decay can then be obtained from E{t) oc which yields 

E{t) (X £P-\t) oc {t - to)~^^ . (8) 

Of course, the result E{t) oc iP~^{t) could also be deduced directly from (^) and the spectral 
model @. For any choice of p, the above spectral decay law is self-similar. In fact, using 
v'^{t)£{t) oc eP-^{t), C{t) oc eP{t), and e'^/^{t) oc (P'^'^/^it), which follow from the previous 
relations, we obtain (||) with 

{k"^ k<1 

The results in (0) and (js]) may be compared with those obtained from numerical solution 
of EDQNM, l{t) oc (t - tof-^^ and E{t) oc (t - to)~^'^^ (see section 10.2 of Lesieur^). From 
this we may infer the value p = 0.55 for EDQNM. However, the above exponent values do not 
depend too sensitively on the precise value of p, which is thus poorly determined by them. For 
example, the value p = 1 implies the relations t{t) oc {t — to)^^^ and E{t) oc (t — tQ)~^/^, which 
are also in close agreement with the numerical EDQNM results. 

The decay laws which follow from this spectral model, for any p, are the same as those 
which are traditionally believed to hold in a self-similar decay for the spectrum with low- 
wavenumber exponent nc = 4 — p and constant coefficient A. This constancy of the coefficient 
A for —1 < n < 4 is called the permanence of the large eddies (PLE)'^. In fact, according to the 
theory based upon self-similarity and PLE, it is usually deduced that 

e{t)cx{t-tor, a=^—, (9) 
n + 3 

and 

Eit)^{t-tor\ 6 = ^^^^ (10) 

n + 6 
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for — 1 < n < 4 (see section 10.2 of Lesieur^, section 7.7 of Prisch^, or Clark and Zemach^). 
We find, if we take n = A— p, that we obtain the same values of a and b as before. There is of 
course no contradiction here, because the coefficient C{t) of the term in (2) depends upon 
time and this changes the decay laws. 

It is no coincidence that the decay laws forn = 4 coincide with those for the standard model 
based upon PLE when n = ric- In fact, we argue, following Gurbatov et al.^ and Eyink and 
Xin^, that there is no self-similar decay at all when Uc < n < 4. We continue to adopt the PLE 
hypothesis. However, we question the additional, implicit assumption in deriving decay laws 
that i{t) is the only length scale in the problem. We consider the consequences of the C{t)k^ 
backtransfer spectrum, with C{t) oc ff'{t) as determined above, and consider the possibility that 
an intermediate spectral range may form which is dominated by the backtransfer. Thus, we 
take as our model spectrum 

Ak'^ kl^{t) < 1 

E{k,t)c^{ c(t)k^ k£4t)>l,ke{t) <1 

e2/3(i)fc-5/3 > I 

Continuity of the spectrum at the juncture = 1 requires that 

ef{t) (X c{t) (X ep{t) 

and thus 

£4t)(x£^{t). 

The intermediate k^ range only survives — and grows — if S> £{t) for large t. Clearly this 
requires that p/{A — n) > 1 or n > 4 — p = Uc- The growth rate of becomes infinitely 
large as n approaches 4 from below, and, in that limit, the k^ region grows to infinite extent. 



For n > 4 — as in the traditional view — the PLE hypothesis breaks down and the decay is 
self-similar, governed by the first model (Q). 

We see that in the range tt-c < n < 4, the decay is not self-similar, as there are two distinctive 
length-scales i{t) and i^{t). The Ak"^ low-wavenumber spectrum does not dominate the C{t)k'^ 
backtransfer spectrum over the whole range ki{t) <C 1, because the latter has an increasing 
coefficient. Instead, the k'^ region is growing in extent and it is this region which matches onto 
the energy range at ki{t) ~ 1. It follows that the decay laws are those determined by the 
backtransfer spectrum @ over the whole range of low-wavenumber exponents greater than nc- 

2 2(5-p) 

£{t) oc {t-to)''~p, E{t) (X {t-to) 7-p n>nc. 

Thus, the decay laws (^), ( p!o| ) hold only for — 1 < n < ric. In our opinion, this is a state of 
affairs far more a priori plausible than the traditionally presented one. In fact, in the traditional 
view there is a monotonic growth in the energy decay exponent b = 2{n + l)/(n + 3) over the 
range — 1 < n < 4, which then suffers a discontinuous drop at n = 4. This seems unphysical. 
In our theory, the exponent b increases as a function of n over the range — 1 < n < ric, but for 
larger n sticks at the value for n = Uc- Although the lowest wavenumber spectrum then does 
satisfy PLE, it does not match onto the energy range and it plays no role in the energetics of 
the decay. In fact, the energy E^,[t) in the lowest wavenumber range where PLE holds is 

Since the total energy scales as E{t) oc for nc = 4 — p, there is a negligible fraction 

of energy in the PLE range asymptotically in time for n > Uc- Hence, at these very long times, 
self-similarity is effectively restored, described again by the first model (Q). 

We presume these facts will also follow from analytical closures, such as EDQNM, since the 
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basic ingredients of our phenomenological theory are already present there. Thus, we expect 
that those closures have no self-similar decay solutions in the range ric < n < A. It is interesting 
that the verifications of self-similarity which have been made do not seem to have included 
values in this range. See for example figure VII-8 in Lesieur^, where an impressive similarity 
collapse is shown, but only forn = 2 and n = 4. We expect that the assumption of a self-similar 
decay in EDQNM-type closures for Wc < w < 4 will lead to a realizability violation, similar to 
what was found for the Kraichnan model^. Of course, positive spectra are guaranteed for 
EDQNM, but only if one actually solves the model and not if one makes hypotheses (such as 
self-similarity) which may be inconsistent with the closure equations themselves! 

A theory analogous to that presented here can be developed for many other turbulent decay 
problems. For example, the case of stationary turbulence with a Richardson eddy-diffusivity 
K{r) ~ r^/^ was studied by Eyink and Xin^, within the Kraichnan model. It was found that 
a decaying scalar in the inertial-convective interval experiences self-similar decay only for low- 
wavenumber scalar spectral exponents which are initially in the ranges — l<n<8/3 and n > 4. 
In this model, the time-dependence of the constant C{t) could be evaluated, as C{t) oc {t — to)^, 
which allowed the determination of nc = 8/3. In the range 8/3 < n < 4, the self-similar decay 
is nonrealizable and is replaced by a two-scale decay of the type described above. This state 
of affairs presumably holds as well for the passive scalar advected by an actual (not synthetic) 
turbulent velocity. In fact, one of us^ has constructed a simple model of the mandoline geometry 
often used experimentally to study decay of temperature fluctuations, and found that there is 
a low wavenumber k'^ spectrum but with decay exponents the same as for n = 8/3. Similar 
results can be derived also for scalars passively advected by turbulence which is itself decaying, 
and for many other such turbulent decay problems. 
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